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Abstract. We reply to comments by P.Markos, L.Schweitzer and M.Weyrauch 
[preceding paper] on our recent paper [J. Phys.: Condens. Matter 63, 13777 (2002)]. 
We demonstrate that our quite different viewpoints stem for the different physical 
assumptions made prior to the choice of the mathematical formalism. The authors of 
the Gomment expect a priori to see a single thermodynamic phase while our approach 
is capable of detecting co-existence of distinct pure phases. The limitations of the 
transfer matrix techniques for the multi-dimensional Anderson localization problem 
are discussed. 
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In our original publication PP a conceptually new view on the two-dimensional 
Anderson localization has been put forward. The main message of our work can be 
summarized in one sentence: the metal-insulator transition (MIT) 'should be interpreted 
as a first-order phase transition' fTl. We believe that a misunderstanding of this crucial 
point (which indeed 'contradicts standard wisdom' has lead to the objections of P. 
Markos et a/.pj. In this Reply we elaborate on this conceptual difference and then 
comment briefly on specific points of the criticism mentioned in [2]. 

One of the cornerstones of the famous scaling theory of localization [Hj and the 
subsequent numerical studies is that for each particular energy and disorder strength 
the system must be only in one of the two phases: with either localized, or extended 
wave functions. The main conclusion of our work does not fit this phenomenological 
framework: we argue not for the existence of a single extended phase, but for the co- 
existence of extended and localized phases for the same energy and disorder strength. 
Which of these two pure phases is realized, depends on particular reahzation of the 
disorder. 

How does this conceptual difference translate into the difference between our results 
and the results of conventional transfer matrix calculations Both approaches start 

with the formal statistical ensemble which includes all realizations of the disorder and 
thus describes both pure phases, extended and locahzed. Averaging a physical quantity 
over this ensemble produces an average between the properties of two pure phases which 
is not capable of characterizing each phase separately. This is the source of failure for 
the standard approach: it considers such heterogeneous averages only and views them 
as characterizing a single pseudo-pure phase. In such a mixture, Lyapunov exponents 
of the pure insulating phase always dominate and this gives a false impression that 
only localized solutions are present. This is how P. Markos et al. pose the question: 
whether the average transfer matrix T^^^ describes the extended or the insulating phase? 
Subsequently, the answer they find is not satisfactory: the assumed pseudo-pure phase 
appears to be localized in 2D (as well as in 3D) for any degree of the disorder. 

Distinctive features of our mathematical formalism (signal theory) follow naturally 
from the requirement that the theory must he able to describe adequately multiple pure 
phases (in case if more than one such phase exist). We have repeatedly emphasized this 
point in our original work by introducing the term 'multiplicity of solutions' (MS) pp. 
The original problem is linear and the MS requirement for an exact solution puts severe 
restrictions on the form of the mathematical formalism. Two points are essential here: 
(i) a precise definition of the phase in analytically exact theories is possible only in the 
thermodynamic limit jH]. That is why our approach assumes an infinite system in the 
lateral direction (M = oo in terms of 0) from the very beginning, (ii) We go beyond 
the description of heterogeneous averages {il'nm'^n'm') (signals in our terminology) by 
considering the fundamental function — the filter H{z). In the region of phase co- 
existence the filter shows MS, corresponding to multiple subsets of propagating signals. 
This is the point where one can start analyzing the properties of distinct pure phases 
and not the mixture of them. 



Reply to Comment 



3 



The transfer matrix approach |2llllini does not posses these crucial features. It starts 
with a quasi-lD system and approaches the thermodynamic hmit only asymptotically. 
We have already emphasized in our work (see Sec. 4.2 in P) that the quasi-lD transfer 
matrices loose the MS property which is the key for a correct description of multiple 
phases. Indeed, P. Markos et al. note [2] that the thermodynamic limit is problematic 
in their approach: 'the limit M oo of the discrete model discussed here bears various 
conceptual and technical difficulties.' Their suggested solution (studying this limit 
numerically) is hardly adequate to overcome these difficulties, while our approach treats 
the system as a truly multi-dimensional one from the starting point. 

It is not a coincidence that we were able to solve the ID case both by the method of 
transfer matrix and the signal theory, while for 2D the former method was abandoned 
[1]. In ID, there is no phase transition and the MS property plays no role, while in 
higher dimensions the two methods are not equivalent. 

Even within the transfer matrix approach, which should give the properties of 
a single pseudo-phase, we do not fully agree with the authors of the Comment 
In the transfer matrix approach (as well as in the signal theory), the fundamental 
Lyapunov exponent is determined by the maximal eigenvalue of the transfer matrix 
(respectively, the maximal root of the equation for filter poles). Therefore the analysis 
of the eigenvalues close to the unit circle j2] is not sufficient to determine the phase 
diagram of the system. 

We can also not agree with the statement (2j concerning the ambiguity of our 
averaging procedure. The latter has been described at length in Sec. 3.2 of our original 
work^, assuming familiarity of the reader with only basic aspects of the signal theory 
and linear algebra. 

Regarding the last critical remark in |2] we note that since our approach does not 
seek to 'detect a metallic phase' but rather to find the multiplicity and the properties 
of pure phases, there is no physical reason to prefer the analysis of the (In I?/'!) over the 
analysis of the second moment (IV'P)- 

Finally, we note that we are not alone in challenging the prevailing view that 
there is no MIT transition for 2D Anderson hamiltonian. The support from finite-size 
scaling studies has been put under the question mark E] • An even more important 
challenge comes from experiment. A MIT has been observed experimentally in 2D 
samples jH], causing significant new research activity. Direct electrostatic probing ITU] 
and photoluminescence spectroscopy ^Tj show a co-existence of localized and metallic 
regions associated with 2D MIT, and new theories are put forward to address this 
issue |l2|il2j. Ref. [THj associates the phase separation with a first-order phase transition 
between a Fermi liquid and a Wigner crystal. At the same time there is a growing 
evidence that the transition is disorder- driven and does not stem from electron-electron 
interactions Thus there is a clear need for a revision of the canonical point of 

view on the localization problem. This revision should touch not only the numerical 
studies of a tight-binding hamiltonian [ZJIH], but mainly the scaling theory of localization 
which stands on phenomenological grounds. Our exact analytic results for a microscopic 
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model show that one does not need to go beyond the original framework of the Anderson 
hamiltonian in order to describe a MIT and phase co-existence in 2D. 
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